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ABSTRACT 

We consider a model in which a collective state couples to a large number of 
background states. The background states can be chosen to have properties which 
are classically characterized as regular or chaotic. We found that the dynamical 
nature of the background system considerably affects some fluctuation properties 
of the strength function. 



1. Introduction 

Over the last twenty years a number of new giant resonances have been found 
in a broad range of the nuclear table. They are coherent particle-hole excitations 
carrying a large fraction of the energy weighted sum rules. In most cases they are 
embedded in a continuum and are damped due to particle escape. At not much 
high excitation energies, however, the damping of giant resonances is dominated by 
spreading, i.e., by coupling to a huge number of background states. It is also known 
that there is in some cases a mechanism which reduces the width of resonances [1] . 

In many theoretical calculations the spreading of a collective particle-hole state 
is accounted for by coupling to more complicated states, e.g., two-particle two-hole 
states. As the number of these states in heavy nuclei becomes enormously large it 
is necessary to truncate the configuration space, or to introduce some kind of sta- 
tistical treatment. For instance, in ref.2, the spreading was discussed by means of 
the random matrix theory. Aside from the reduction mechanism mentioned above, 
the coupling matrix element to each of the background states is expected to be- 
have randomly, reflecting a random nature of the latter. In fact, the single-particle 
width distribution of neutron resonances obeys the Porter- Thomas distribution 
consistent with the random matrix theory [3]. We may note that level statistics 
of excited states have been investigated in realistic shell model calculations[4,5]. 
Not much is known, however, if or what properties of a strength distribution of 
collective state are related to the dynamical properties of the background states. 

In the present note we consider a model which imitates a coupling of a collective 
state to background states, the latter being chosen to have chaotic or regular 
dynamical properties. In this manner we investigate wheter there is a characteristic 
signature in the strength function which reflects a nature of the background system. 

2. Model 

In the study of giant resonances, the Hamiltonian is normally diagonalized 
within the truncated space of one-particle one-hole states plus background states. 
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mainly two-particle two-hole states. In the present model we replace the back- 
ground states with eigenstates of a system whose classical counterpart is well stud- 
ied, while a collective mode is represented by a single boson excitation independent 
of the background system. The structure of the collective strength function de- 
pends both on the coupling Hamiltonian of the collective mode to the background 
and also on the dynamics of the background system itself. We adopt a simple 
ansatz for the coupling in order to focus on the latter effect. We start with the 
Hamiltonian 



H = i?coll + -f^bg + Koupb 



(1) 
(2) 



Here //coll represents a Hamiltonian for the collective excitation, where B'^ and B 
represent boson creation and annihilation operators. We write the boson vacuum 
as |0). As the Hamiltonian i/^g of the background system we choose that of a 
two-dimensional anharmonic oscillator characterized by a single parameter k. The 
classical counterpart of this Hamiltonian describes a transition from an integrable 
to a chaotic dynamical system[6]. The quantum spectra and the wave function 
characteristics follow the same trend [7, 8]: At A; = the system is separable while 
at = 0.6, for instance, the nearest neighbor spacing distribution of the quantum 
levels shows the Wigner distribution which is typical for chaotic systems. We 
denote the eigenstate of i/^g as ( \n = 0)) for the ground state). In the 
diagonalization of -f^bg; we take as the basis states of the background system the 
eigenstates of an uncoupled harmonic oscillator whose frequency fl is determined so 
as to optimize the diagonalization of i^bg [7,8]. They are denoted by \N) (\N — 0) 
for the ground state) where N stands for a pair of integers, i.e., numbers of oscillator 
quanta in the x— and the y— directions. 

The interaction Koupl represents the coupling between the collective state and 
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the background states. We take a simple ansatz for the couphng, 

V,,^p, = xJ2iB^\0){N\+h.c.) (3) 

characterized by the strength parameter x- Note that the state \N) is different for 
each value of k because of the difference in the Q value. Due to the simple form 
of eq.(3), regular or random behavior of the coupling matrix elements depends 
entirely on the dynamics of the background system which can be controlled by the 
parameter k. The assumption (3) may not be realistic, as the coupling strength of 
the actual nuclear system would decrease for complicated states (i.e., large N). We 
expect, though, eq.(3) is sufficient for the present purpose of studying a qualitative 
difference coming from the dynamical structure of the background system. 

We divide the whole space into four parts: 

1. The ground state of i^coll + ^bg: |0; 0) = |0)|0)). 

2. The one collective boson state: |1; 0) = B^\0)\0)). 

3. The background states: |0;n) = |0)|n(7^ 0))). 

4. Other states. 

In diagonalizing the Hamiltonian (1), we neglect the coupling between space 1 and 
space 2, and also omit the space 4. These approximations produce negligible effect. 

The eigenstates of i?bg are classified into several symmetry classes with no 
coupling among them [7]. In the present calculation we consider only those states 
which belong to the class symmetric in the x—,y— and the diagonal (a; = y)- 
directions. First i/bg is diagonalizcd within a large space (the number A^max of basis 
states is 5776) and then the lowest 800 states are included in the diagonalization of 
the total Hamiltonian. As for the values of k, we consider three typical cases, i.e., 
k —0.0 for an integrable background system, 0.2 for a partially irregular system, 
and 0.6 for an almost chaotic system. The value of e has been fixed to 220, so 
that the collective state is located in the middle of the background 800 states. 
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and thus a large number of background states can be found in the neighborhood. 
For a fair comparison of the role of integrable versus chaotic background system 
one might better adjust e value for each A; to give a similar local level density. 
Although the asymmetry in the level density of the background states around e 
does affect the third moment of the strength function as discussed later, its effect 
on the fluctuation properties of the strength distribution is expected to be small. 
The coupling strength is mostly fixed to x = 1-0 in the calculation below. 

3. Distribution of Strength 

Collective strengths are measured with respect to the operator O = + B, 

i.e., 

S{E) =J2^iE- {Ea - Eg.s.))\tot{a\ 6 \g.s.),^, \\ (4) 

a 

where |q;)^q^ denotes an eigenstate of and E^ the corresponding eigenvalue. 

Before discussing the strength function S{E) in detail, we first examine the 
distribution of the coupling matrix elements Vn = (1; 0| Koupl |0; ^t-)- According to 
our choice of the interaction it is expected that the coupling matrix elements Vn 
would behave regularly for k small and randomly for k in the chaotic regime. This 
is indeed so as seen in Fig.l where the distribution of the coupling matrix elements 
is shown. 



Figure 1 



It is seen that the matrix element values for k =0.0 are concentrated at ±1, while 
for k — 0.6 they show almost the gaussian distribution centered at zero and with 
the width of around 1.1. In fact, for a variable Vn composed of a large number of 
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Figure 2 



random elements the distribution would be a gaussian, 



-\/27ra; 



1 



'Tl 



exp( 




(5) 



with a-n — 1 from the normalization of Vn- 

Figure 2 shows a strengh function S{E). In spite of the difference in the 
coupling matrix elements as seen above, the shape of S{E) looks rather similar for 
k=0.0 and 0.6. For the sake of quantitative discussion of the gross structure, let us 
consider the cumulant {E"^)^ of the strength function. Several low order cumulants 
are given by 



{E)^^{E), {E^)^^{E^)-{E)\ {E')^^{E')-3{E^){E)+2{E)\ (6) 



where (£'" ) is the n-th energy-moment of the strength function defined as, 



We show in Table 1 the calculated cumulants up to the third order at three k values 
and for % = 1.0. Let us consider the second cumulant which can be written as. 




(7) 




n(#0) 



(8) 




iV,iV'(#0) 



Table 1 
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Eq.(8) shows that, if we do not truncate the number of |n)), we would obtain almost 
the same values of {E^}^ determined by A^max for any k [9]. Indeed, the actual 
values of (-E^)^, as listed in Table 1, are all close to 800 for the three values of k. 
The third cumulant can be written as, 

= E K - '>n, (9) 
n(^0) 

where Un denotes the eigenvalue of i?bg- From Table 1, we find that the value 
of (E^)^ decreases as k increases. This trend is mostly due to the difference in 
the level densities around e as mentioned before: If we artificially use the same 
sequence of cun for the three k values we obtain almost the same values for the 
third cumulant. The apparent similarity of the gross structure of the strength 
function may be understood in this manner. We mention that the shape of the 
strengh function is not much different for different values of x, except that the 
width of the distribution is scaled accordingly. 



Figure 3 



The similarity of the strengh function is only superficial, however. This is 
seen in the strength distribution P{S) as given in the upper part of Fig. 3. The 
smooth curve shows a Porter-Thomas distribution which is expected to hold in 
the chaotic system as given in the random matrix theory [3,4]. Difference due 
to the dynamics of the background system can be more clearly seen in the lower 
part of Fig. 3, where the distribution of the amplitude \/S corrected for the energy 
dependent factor ((i? — e)^ + is shown. The latter has been introduced to 
remove approximately the energy denominator contribution: If we assume constant 
coupling matrix elements Vn = Vf, and an equal level distance D^-, the strength 
function will be given by S{E) ~ (r/27r)/{(£; - e)^ + (r/2)2}, where F = ^t^vI/D^ 
[9] . This may be contrasted to the present calculation where the mean square value 
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of Vn gives ~ X^(see Table 1) and the mean level distance is D ~ 0.5. Thus 
we choose F = Attx in Fig.3. The result shows that at A; = 0.6 the distribution 
follows a gaussian, while at k = 0.0 it is peaked at unity. Notice that the former 
distribution is generic, while the latter, corresponding to the regular background 
system, depends on the specific structure of the coupling. 

The strength distribution is one of the signatures which characterize the un- 
derlying dynamics of the system. It reflects only a part of the structure in the 
strength function of Fig. 2. For instance, the structure of the eigenstatcs of i/bg 
as a function of energy changes regularly at small k values, while at A; = 0.6 it 
strongly fluctuates from state to state[8]. A similar behavior may be expected to 
occur in the present coupled system. For this purpose we perform a moment anal- 
ysis similar to the one applied for a multifractal system [10]. This analysis takes 
into account some features of the energy-strength correlation, and therefore can be 
another characteristic measure of the strength function independent of the distri- 
bution P{S). Here the whole energy interval AE is divided into L segments each 
having a width SE = AE/L. (We take AE = -Egooth — Eg,s. in the present case.) 
The strength in each segment is summed up to give the strength Pj{j = 1, . . . ,L) 
for the j— th segment, with the normalization condition Y2j — 1- The m— th 
moment at the scale SE is defined by 

L 

Mm{5E) = J2Pr- (10) 
We then study its behavior as we refine the scale, e.g., as L = 2 —>■ 2^ —>■ 2^ etc. 



Figure 4 



Figure 4 shows the dependence of the moments Mm for m = 2 to 5 on the 
scale SE at k — 0.0 and 0.6. We also show, for the sake of comparison, the one 
for the ideal case of equidistant eigenvalues with equal strength. Since we have 
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a discrete spectrum, the moment Mm for small 5E, i.e., for a large number of 
segments eventually reaches a fixed value. It is seen that the slope becomes almost 
constant around dE c^2. We now consider the fractal dimension Dm defined by 

Dm = lim — , Bm{oE) = — . (11) 

5E^Q m — 1 logSE 

In practice, the expression for Bm{SE) is replaced with the ratio of the difference 
of log Mm{SE) to that of logSE in the appropriate interval of 5E. The quantity 
Dm reflects a state-to-state fluctuation of the strength function. In Table 2 we 
show the calculated values of Dm for m = 2 to 5. 

Table 2 



The adopted interval of 6E is indicated in Fig. 4. The result shows that the strength 
function at k =0.6 has a smaller Dm value than that for k =0.0 and for the ideal 
case. We also see from Fig.4 that the dependence of Mm on SE is smoother 
for k — 0.6 than the other cases. These results may indicate that the scaled 
moments (10) provide another characteristic signature which reflects the underlying 
dynamics of the system. 



4. Summary 

In summary, we constructed a model which simulates a collective state cou- 
pled to a background. We took a two-dimensional anharmonic oscillator as the 
background system which exhibits a regular or a chaotic spectrum depending on 
the parameter k. We found that the difference in the dynamics of the background 
system causes a characteristic difference in the strength distribution. We also 
suggested that the scaled moment analysis might be useful in the study of an 
energy-strength correlation, although it requires a further study in order to make 
clear whether the conclusion survives in a more generic context. In the actual 
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study of the strength function of a physical system, e.g., of the nuclear giant reso- 
nances, the finite experimental resolution as well as the effect of a continuum put 
a restriction on the applicability of the method presented here. Other possible 
methods, such as the autocorrelation function and the Fourier transform, are now 
under investigation [11]. 

References 

[1] G.F.Bertsch, P.F.Bortignon and R.A.Brogha, Rev.Mod.Phys.55(1983)287. 

[2] T.Guhr and H.A.Weidenmiiller, Ann.Phys. 193(1989)472. 

[3] C.E.Porter and R.G.Thomas, Phys.Rev.l04(1956)483. 

[4] T.A.Brody et al., Rev.Mod.Phys. 53(1981)385. 

[5] M.S.Bae et al., Phys.Rev.Lett.69(1992)2349. 

[6] H.-D. Meyer, J.Chem.Phys.84(1986)3147. 

[7] Th.Zimmermann et al, Phys.Rev.A33(1986)4334; 

O.Bohigas, S.Tomsovic and D.UUmo, Phys.Rep.223(1993)43. 

[8] H.Aiba and T.Suzuki, RCNP 1991 Annual Report and to be pubhshed. 

[9] A.Bohr and B.R.Mottelson, Nuclear StructureVol.l (Benjamin, 1969) Chap.2D.| 
[10] J.L.McCauley, Phys.Rep.l89(1990)225; 

see also, S.N.Evangelou and E.N.Economou, Phys.Lett.A151(1990)345. 
[11] H.Aiba, S.Mizutori and T.Suzuki, in preparation. 



10 



TABLES 



Table 1. Cumulants of the strength function from the first to the third order for % = 
1.0. 



k 








0.0 


220.0 


842.5 


28516.5 


0.2 


220.0 


779.5 


15784.9 


0.6 


220.0 


965.9 


-1942.3 



Table 2. Fractal dimension D^n (m = 2 to 5) at A; = 0.0 and 0.6 for x = l-O- Also 
those for the ideal case of equidistant energy eigenvalues with equal strength 
are listed for comparison. 



k 


D2 


^3 


D4 


D5 


0.0 


0.73 


0.71 


0.70 


0.69 


0.6 


0.64 


0.56 


0.51 


0.47 


Ideal case 


0.98 


0.98 


0.97 


0.97 
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FIGURE CAPTIONS 



Fig.l Distribution of the coupling matrix elements Vn for x = 1-0 at three values 
of k. The smooth curves show a normalized gaussian distribution having the 
same width as that of S{E) at each k value. 

Fig.2 Strength function S{E) for x = 1-0 at k =0.0 and 0.6. 

Fig.3 Upper part: Strength distribution P{S) for x = 1-0 at three values of k. The 
smooth curves show the Porter-Thomas distribution. Lower part: Distribu- 
tion P{^/S) where S ^ S{{E - + {T/2f). T is fixed to 4n. The smooth 
curves show the gaussian distribution. 

Fig A Mm{6E) versus 5E for x = 1-0 at k =0.0 and 0.6 as well as for the ideal case 
of equidistant energies with equal strengths. The lines correspond to m =2 
to 5 from the upper to the lower ones. The arrows indicate the interval of 
SE where the fractal dimensions are evaluated. 
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